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Abstract 

The charmonium system is studied in a Salpeter model with a vector plus 
scalar potential. We use a kinematical formalism based on the one devel- 
oped by Suttorp and present general eigenvalue equations and expressions 
for decay observables in an onium system for such a potential both in the 
Feynman and Coulomb gauges. Special attention is paid to the problem with 
renormalization of the lepton pair decays, and we argue that they must be 
defined relative to one of the experimental decay widths because renormal- 
ization of the vertex function is not possible. The parameters of the model 
are determined by a fit to the mass spectrum and the lepton pair decay rates. 
Two gamma decays and El and Ml transitions are then calculated and found 
to be well accounted for. No significant differences in the results in Feynman 
or Coulomb gauge are found. A comparison is made, regarding the electro- 
magnetic transitions, between the full and reduced Salpeter equation. A large 
difference is found showing the importance of using the full Salpeter equation. 
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I. INTRODUCTION 



The Bethe-Salpeter formalism has recently been used in various ways in more careful 
investigations of quarkonium physics. For example Munczek and Jain have used the 
full Bethe-Salpeter equation together with the Schwinger- Dyson equation to study the mass 
spectrum and lepton pair decays of light and medium-heavy mesons, Resag et al. 
have written a series of papers about the Salpeter equation where they have studied mass 
spectra and various decay observables for both light and heavy mesons, Kugo et al. have 
calculated the Isgur-Wise function and Harada and Yoshida have calculated the QCD 5* 
parameter and developed a method to solve the full Bethe-Salpeter equation 0. The 
stability of the Salpeter equation for various Lorentz structures of the confining potential 
has been studied by Olsson et al. [|TU| and Parramore et al. |]TT],|T^. Olsson et al. have also 



studied the validity of the reduced Salpeter equation |T3|. 

In this paper we undertake an investigation of the radiative decays of the lowest char- 
monium levels including both Ml and El transitions in our study as well as the two-gamma 
decays of the corresponding pseudoscalar, scalar and tensor states. 

Like many other authors [p|, p!^p!5| we use the instantaneous approximation of the Bethe- 
Salpeter equation - the Salpeter equation - which has a well defined physical interpretation. 
Our kinematical analysis uses a formalism based on the one developed by Suttorp [jl| . (An- 



other formalism in a similar spirit has been developed by Lagae [|16|.) This leads to a set 



of coupled equations for the amplitudes which are solved numerically. In this paper we 
present a general set of eigenvalue equations and decay amplitudes for lepton pair decays, 
two gamma decays and El and Ml transitions, which can be used in any model with a 
vector plus scalar potential using the formalism of Suttorp. 

Our potential consists of one short-range vector part describing the one-gluon-exchange 
interaction and one scalar long-range confining part. Investigations of the stability of the 
Salpeter equation for confining potentials with various Lorentz structures have shown that for 
light quarks there are no stable solutions for a scalar confinement [p!0|-p!^ . This problem does 
not appear for heavy quarks like the charm quark. Furthermore a Lorentz time component 
vector confinement, which has been shown to be the most stable one, gives a contribution to 
the spin-orbit coupling with the wrong sign. We have therefore chosen a scalar confinement 
potential. 

In quantum mechanics the spectrum itself does not uniquely determine the potential. 
We also need as input some other information concerning the localization of the system, 
such as the value of the wave function close to the origin. As input we therefore take the 
energy levels of the charmonium spectrum as well as the decay rates of the vector 1 states 
into lepton pairs. 

The instantaneous approximation leads, strictly speaking, to some problems with renor- 
malization, since the Ward-Takahashi identities are then not relevant for the renormalization 
of the vertex function in the lepton pair decay. Furthermore the kernel is confining and con- 
tains a dimensionful coupling constant for the confining part. This leads in the final analysis 
to a renormalization procedure where we use the decay width for J/\E' — e^e~ to define 
the renormalization constant. Thus, we calculate the ratios of the decay rates of the vector 
states into lepton pairs, rather than the rates themselves. 

We adopt the formalism of Mandelstam il^ to calculate the radiative decay amplitudes. 
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The gauge invariance of this method in the instantaneous approximation is not obvious, but 
in Sec. ^ we show that the decay amphtudes are indeed gauge invariant. 

The paper on heavy quarkonia by Resag and Miinz overlaps to some extent with our 
paper. They have also used a Salpeter model to calculate mass spectra and decay observables 
for charmonium. However, they have used a somewhat different potential. In particular, 
they have introduced a cut-off in the potential at the origin in order to renormalize the lepton 
pair decay widths. In this paper we discuss the difficulties of renormalization and conclude 
that an experimental decay width has to be used to properly define the renormalization. A 
cut-off also necessitates some extra parameters in the model, and they obviously affect the 
results. In particular, the lepton pair decay widths depend on the wave function's value at 
the origin, and therefore one can obtain any value for those by varying the cut-off parameters. 

In addition, we present general equations in a formalism different from and investigate 
the importance of non-dominant terms in the wave functions. 

The outline of our paper is as follows. In Section y we briefly discuss the Bethe-Salpeter 
equation to define our variables. In Section |T| we describe the potential used in the model. 
We use a short-range Coulomb potential plus a long-range confining linear potential. In 
Section |^ we give the expression for lepton pair decay and discuss the problem with renor- 
malization. In Section |V| we give the expressions for the 77 decays and in Section [V| we give 
the expressions for the Ml and El transitions. There we also discuss the gauge invariance 
in the instantaneous approximation. In Section |VII| we present the numerical results and 
finally in Section |^ we give a summary of our results and conclusions. 



II. THE SALPETER EQUATION 

The wave function x(g) for bound states of a fermion-antifermion pair satisfies the spinor 
Bethe-Salpeter equation 

+ f - ^)x{qM -^-m)=^j:J ^,K\p, g, P)Px(p)^^ (1) 

Here, P is the momentum of the bound state and q is the relative momentum of the con- 
stituent fermions. The interaction kernel is characterized by the matrices F*, where F"^ = 1, 
YV ^ F^ = F^ = 7V, = 7^- 

We now make the instantaneous approximation, so that the kernels = K^{p — q) 
in the rest frame. This reduces the Bethe-Salpeter equation to the Salpeter equation. It 
may seem like Lorentz invariance is lost by making this approximation. However, although 
manifest Lorentz invariance is lost, the Salpeter equation can be formulated covariantly by 
letting depend on = g — jj^P and p± = p — f^-P in a general frame. In the rest frame 
q± reduces to (0, q). For details, see for example Chang and Chen [p^ . 

It is most convenient to calculate the mass spectrum and wave functions in the rest frame 
of the bound states, and afterwards Lorentz transform the wave functions when needed. We 
will therefore perform the calculations in the rest frame for the remainder of the paper. 

The propagators can be written 
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m \{qO - E{q) +ie) {q^ + E{q) - le) J ' 







^ + ^ ], (2) 



where 



Aj(q) = ^^(i?(q)+/i(±«-q + ^7°)) /^=± (3) 

are the projection operators for the positive and negative energy states respectively. 

With the help of these projection operators we can write the Salpeter equation as a set 
of four coupled equations, 

(q)x(g)A, (q) = (m/2 + q^ - ^lE + fiie){-M/2 + q^ - uE + me) ^ 

E / (^^^(P - q)A^(q)/rx(p)rVA;(q). (4) 

The right hand side can be integrated over g° by a simple contour integration. Introducing 
the reduced wave function 

^(q) = / ^X{q) (5) 

we find that 

A+(q)*(q)A:(q) = - ^ E / (^^^(P " q)A+(q)7°rvl/(p)rVA:(q), (6) 

(q)M/(q)A;(q) = ^^^^ E / J^s^'^P ' q)Al (q)7°r^(p)rVA;(q), (7) 

A+(q)v^(q)A;(q)=0, (8) 
Al(q)^(q)A:(q)=0. (9) 

Combining these equations with (^), we can write the wave function as 

xiq) = 

i{M - 2E) 



(M/2 + gO - E + ie){-M/2 + q^ + E - ie 



-A+(q)^(q)A:(q) 



+ (M/o^ 0^ F '^^u'^fJ/o^ ^7— -TAl(q)^(q)A;(q) 

(M/2 + q'^ + E — ie)[—M/2 + q" — E + le) 



+ 7 -, n 77 ; ^ 7 / ^-^^*A+(q)7"r^(p)rVA7(q) 

(M/2 + g0-E + ie)(-M/2 + g0-E + ze) y (27r)3 ' ' 

+ , n , ^ , n , ^ ■^ / T^i^^A^ (q)7°rv]/(p)rVA:(q). (10) 



(M/2 + gO + E - ie)(-M/2 + gO + E - ze) 7 (27r)3 
Thus, finding the wave function x(g) in this formalism reduces to finding \[^(q). 
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Combining the equations (||)-(H) we get an equation for \l/(q), 
(-^ a • q - m7'^)w(q) + W(q)(— - a ■ q + niY) 

= -E / (^^Xp-q)(A:(q)7°r^(p)rVA:(q) -Ai(q)7°rvi/(p)rVA;(q)). (ii) 

For each state with a certain angular momentum, parity, and charge parity, we expand 
the function \l/(q) in the Dirac algebra by imposing the correct transformation properties. 
We also demand that \I'(q) fulfills the conditions For example the expansion of the 

O""*" state is 

vl/(q) = M/p(g)7^ + ^> j,{q) (7V + 7^^) ■ (12) 

The scalar functions ^ a and \l/p depend only on g = |q|. Similar expansions for other states 
are given in Appendix 0. 

Integral equations for the scalar functions can be obtained from (p!TD by expanding 
both sides in the Dirac algebra and comparing the coefficients. This is discussed further in 
Appendix 0, where all integral equations are given. 



III. THE POTENTIAL 



We have chosen a potential which consists of one short-range part describing the one- 
gluon-exchange interaction and one long-range confining part. The calculations have been 
done both in the Feynman gauge and in the Coulomb gauge, in an attempt to study the 
gauge dependence in the instantaneous approximation. The respective potentials are 

irFeyn(r) = {r)^\^h^^^'^ + {v) + (13) 
KUr) = ir^(r) (70(1)70(2) - 1^(1)^(2) + 1(^(1) . x)(^(2) . x)) + K'{r) + U, (14) 

where ?7 is a constant and 

ir^(r) = ^, (15) 
r 

K^(r) = Xr. (16) 
In momentum space the vector part in the Coulomb gauge becomes 

(7?V^^) + (y^)-(p-q))(7^^)-(p-q))) (17) 



(p - q)^ 



The color factor 4/3 is included in ag, c^s = 4/3as. We let the coupling be constant, 
since a running coupling would not make any significant difference for heavy mesons like 
charmonium. 

There are four parameters in the model: ds, A, U, and the quark mass m. These will 
be determined by a least square fit to the mass spectrum and the ratios of the lepton pair 
decays. 
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Olsson et al. |T0| have recently investigated the stabihty of the Salpeter equation for 
various Lorentz structures of the confining kernel. They have come to the conclusion that for 
light and medium-heavy quarks numerically stable solutions exist only for a time component 
vector potential. The same result has also been found by Parramore and Piekarewicz |TT|JT^ 
from a similar analysis. 

However, we have found no such problems for heavy quarks like in the charmonium 
system. This has also been confirmed in Ref. [Q. Furthermore, a time component vector 
confinement gives a contribution to the spin-orbit coupling with the wrong sign, which speaks 
in favor of a scalar confinement. 

We want to be able to perform the angular integrations in the Salpeter equation analyti- 
cally. We therefore need the expansions of the potential in spherical harmonics in momentum 
space as 

K\v> - q) = E FiiP^ q)Y;^{v)Yim{^). ^ = V^S. (18) 
The expansion of the Coulomb potential is 

where Qi are the Legendre functions of second kind. Thus 

fJM^'-^qA'^]. (20) 

pq \ 2pq ) 

In the Coulomb gauge potential ( p!7D there is another angular dependence in the second 
term. However, we later choose to give this dependence in terms of the functions ■ 

The Fourier transform of the linear potential can only be defined in a distributional 
sense. This causes no problems as in momentum space the value of the potential at each 
point is not important, but only its action on a wave function. We define K^[r) as the limit 

KS[r) = \im^{^^\. (21) 

This gives 

ir^(p-q)=lini — 



e->o dt^ v(p ~ q)^ + 

= !!53 E gr^Q' [ ^„ ] yL(p)y.M- (22) 

It can be shown [|19|, that in a distributional sense, the following relation holds: 

'p^ + + e"' 



52 / 




-j fip) 

'^'^ ^ (p+.)'(p-.)- H ' ^''^ 
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where PV means the principal value. Thus 



IV. LEPTON PAIR DECAYS 

We now continue with a discussion of the decay observables. A general prescription for 
calculating matrix elements of the electromagnetic current operator between bound states 



has been given by Mandelstam [|T^. We begin with the lepton pair decays, which we will 
use together with the mass spectrum, in order to optimize the four parameters as, A, U, 
and m in the model. The probability for a vector particle decaying into two light leptons is 

r = ^b%P, (25) 

where is the vector quark current, Cq is the charge of the quark, and M is the mass of the 
vector particle. To lowest order in the electromagnetic coupling the current is, in the rest 
frame of the decaying particle. 

The corresponding Feynman diagram is given in Fig. |I|. 

Taking the trace we find that j° = (whence the current is conserved) and 

j = / ^ (^-^^^ s{q)YUk) + 4vI/v(g)YS(q)^ . (27) 



We have here used the expansion (|A2|) of the 1 state given in Appendix ^ 



Performing the angular integration we get 

dq 



IJ 



2 



Smq^^>s{q)+Sq^^^v{q)U ■ (28) 



(2vr) 

The expression for the leptonic decay is divergent and has to be renormalized |^. This 



can in principle be done in several ways. However, here we must take into account the 
fact that the theory is not fully defined. This means that with an interaction kernel which 
is instantaneous and defined by the ladder expansion only, the Ward-Takahashi identities 
are not appropriate. The divergence of the wavefunction at the origin can therefore only 
be related to the (experimental) value of the vertex function at some suitable momentum 
transfer. For non-confining potentials, such as the Coulomb potential, there is a region of 
production of asymptotically free fermions above threshold. The production cross section 
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can then be calculated with the kernel and the result be taken as a measure of the vertex 



function. This is more or less the choice in Ref. In our case the potential is confining 
and consists of two parts, one of which has a dimensionful coupling constant A. This means 
that the renormalization has to be done in the deep inelastic region, where the quarks are 
asymptotically free and the linear part of the potential is negligible. This presents several 
problems since the lack of manifest covariance makes the renormalization constant depend 
upon both the cut-off and the momentum. 

The problem appears as follows. We would like to define a renormalization constant Zi, 
such that 

.A 

lim Zi(A) / rfq7^^(q) = finite. (29) 

A^oo J 

Asymptotically \E'(q) ~ q^^^^ , and thus ~ A'^ for large A, where (3 is determined from 
the dynamics |^ . 



The renormalization constant, Zi, is defined through the (analytically continued) vertex 
function by 

u{ci)V^u{ci) = Z^\q)uici)ru{ci) (30) 
and the vertex function satisfies 



r(q) = 7 + ^/ TT^^'^lp - q)7^;^ r(p 



7m 



+ M — ^(p);^ — • (31) 

J {27cj^ p — m p — m 

The renormalization constant Zi has the correct asymptotic behavior to cancel the A- 
dependence of the current. However, due to the instantaneous approximation, Zi depends on 
q. With only a Coulomb potential, Zi has the correct A-dependence to cancel the divergence 
of as A ^ oo for any momentum q. However, with a confining part added to the potential 
this would only be true in the region q ^ oo, where the linear potential is negligible. 
Unfortunately, Zi ~ q^ when g — » oo, and therefore the theory is not renormalizable in this 
way. 

The electromagnetic theory can only be used in lowest order due to the non-existence 
of Ward-Takahashi identities in the instantaneous approximation since the theory is not 
fully defined. Nevertheless there is a loop contribution in the current which has to 
be renormalized. As the vertex function cannot be renormalized, we therefore adopt the 
philosophy that, for confining potentials, the most natural experimental datum at which to 
fix the normalization is the resonance formation itself. The renormalization should therefore 
be effected by using, say, the decay width for the lowest resonance as input. The other decay 
widths are then calculated relative to this. Thus we will optimize only the ratios of the decay 
constants for successive nS states as predictable by the model. These ratios are independent 
of the cut-off A for A > 30 GeV to within 1%. 

Another way to avoid the problem with infinities would be to redefine the potential to 
take a finite value for r = 0. This would make the current (^) finite. However, the theory 
must still be renormalized although only a finite renormalization is necessary in this case. 



8 



Therefore such an approach would still face the same problem as explained above, as a new 
dimensionful constant must be introduced in the cut-off of the potential. Our opinion is 
that our approach is much cleaner and conceptually more transparent. 



V. 77 DECAYS 

The calculation of the 77 decays is most easily performed in the rest frame of the bound 
state, which decays into two photons with momenta ki = (M/2, k) and k2 = (M/2, — k). To 
lowest order the matrix element is given by 

The corresponding Feynman diagrams are given in Fig. 0. The decay probability is given by 

2(2J+ 1)M^ ' 1 2 M'^l ) \ J 

where is the charge of the quark, and M is the meson mass and J its angular momentum. 

The amplitude eie2M^^ can be expressed in the scalar amplitudes for the different de- 
caying states. This is discussed in Appendix ^. 

VI. Ml AND El TRANSITIONS 

In the rest frame of the initial particle the matrix element for the Ml and El transitions 
is, in the ladder approximation, given by 

= Tr y - k/2, Pf)YUq, M,){-^^' + i - m) 

+Tr J |^X/(g + ^/2,P/)(^7° + ^-m)x.(g,M.)7^ (34) 

Here XiiQi^i) is the wave function of the initial particle in its rest frame, and i 
k/2,Pf) is the wave function of the final particle in a frame where it moves with momentum 
Pf. This wave function is found by Lorentz transforming the one in the rest frame, see 
Appendix p 1| . The Feynman diagrams corresponding to are given in Fig. ^. 

The decay probability is given by 

^- ' El^/^^l' (35) 



2M2(2J+ 1) 

where is the charge of the quark, k the photon momentum, Mj the mass of the initial 
meson and J its angular momentum. 
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To calculate the matrix element we use the expansions of the wave functions given in 
Appendix 0. We first perform the g° integration. Since the Lorentz transformed state 
x{q — k/2,Pf) has a very complicated analytic structure we cannot make a simple contour 
integration. We have to use the principal value technique. The real part is given by a sum 
of residues 

where is a zero of g{p^)- The imaginary part gives an unphysical non-Hermitian contri- 
bution to the decay width. It should vanish in any consistent covariant theory, but in the 
Salpeter model is has been verified that it does not vanish |0, because of the lack of manifest 
covariance. However, for the charmonium system this unphysical contribution is negligible, 
and we therefore neglect it. 

The integration over the azimuthal angle ip can also be done analytically. The rest of the 
calculation is made numerically. Trying to take the trace analytically creates enormously 
complicated expressions which are hardly of any use. 



A. Gauge invariance of the matrix element 

In a Bethe-Salpeter model without the instantaneous approximation where we have an 
even kernel of the form K{p — q), the gauge invariance of the matrix element can be shown 
with the help of the Bethe-Salpeter equation. Let the initial meson have momentum Pj and 
the final meson have momentum Py = Pj — k, where k is the photon momentum. Then 



k,M^ = J ^,x{q-k/2,Pj)Mq.Pi){-^ + i-m) 



(2vr) 

= / ^,x{q-k/2,Pf){^ + i-m)x{q.Pi){-'^ + i-m) 

" / (If^ (^^^^ ~ ^^^^^^ ~ = ° ^^^^ 

provided the order of integrations can be interchanged. This proves gauge invariance. 

However, if the kernel depends on the momentum of the bound state, as it does in the 
instantaneous approximation, the terms in the last line do not cancel. In the first term K 
depends on Pj, and in the second term on P/. 
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Nevertheless, in the instantaneous approximation we have verified that k ■ M = and 
M° = in the rest frame of the decaying particle, by performing the angular integration of 
the decay amplitude. Thus gauge invariance is fulfilled. 

VII. NUMERICAL RESULTS 

The eigenvalue equations given in Appendix P were solved numerically by expanding the 
wave function in cubic Hermite splines |^. However, first the integration region < p < oo 
was transformed into the finite interval — l<a;<lbya; = (p — po) /{p + Po), where po was 
chosen in the physically relevant region. The integration and the following matrix eigenvalue 
problem for the expansion coefficients were solved with standard methods. We have tested 
the stability of the solutions by varying the number of basis functions and the parameter 
Pq, and found that the eigenvalues do not change. Problems, like the ones reported in Refs. 
||T0HT2| for lighter quarks, have not been seen. 

The parameters a^, m, A, and U of the model are determined by making a least square 
fit to the mass spectrum and the ratios of the lepton pair decay widths J/\l//\l/(25') and 
^(25')/^'(35'). As discussed in Sec. 0we renormalize the decay widths to the experimental 
value of the decay of J/^. We therefore make the fit to the ratios between the widths rather 
than the widths themselves. We minimize the function 

2 ^ (M,(theory) - M,(exp))^ , ^ (F, (theory) /F, (theory) - r,(exp)/r,(exp))^ 
X =2^ -2 + 2^ -2 • l-^^J 

The first sum goes over all the charmonium masses in Table | except t]c{2S) and h^i which 
are not used in the fit because they still need experimental confirmation. The second sum is 
over the two ratios of lepton pair decay widths r(J/^)/r(^(25)) and r(^(25))/r(^(35)). 
We have not used the decay widths for the higher states because their wave functions are 
not good enough to be used for calculations of decay widths. The experimental errors for 
the ratios, ai/j, are numerically aj/<i,/2s = 0.30 and <J2s/3S = 0.64 respectively. 

As many of the masses are very accurately measured, using the true experimental error 
for the masses in this function would make the masses totally dominate over the decay width. 
We therefore use a minimum value of (Xj when the experimental error is small. Two fits for 
both gauges have been made, one where we have used cTj = 10 MeV for the masses when the 
experimental error is smaller than 10 MeV, and one where we have used cxj = 50 MeV for the 
masses when the experimental error is smaller than 50 MeV. We use the true experimental 
errors for ai/j in both fits. In the second fit we thus give the lepton decays a higher weight, 
though in both cases the masses get a higher weight than the decay widths as the errors 
of the decay widths are rather large. The numerical results of the four fits can be found in 
Tables | and In these tables F and C stand for Feynman and Coulomb gauge respectively, 
and 10 and 50 MeV stand for the minimum value of Uj. 

Making four fits gives us a chance to see how sensitive the different decay observables are 
to changes of the model parameters. The parameter values in the four different cases can be 
found in Table |IT1[ We see that when we change ai from 10 MeV to 50 MeV the parameter 
which changes most is the constant U. One would perhaps expect that as would be the 
parameter most sensitive to such a change in cTj, as the lepton pair decays should be very 
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sensitive to dig because they depend on the value of the wavefunction at the origin. However, 
we calculate the ratios between the decay widths rather than the widths themselves, and 
these ratios are not very sensitive to the value of a^. This can be understood by noting that 
in a non-relativistic model with a Coulomb potential the ratios of the wave functions' value 
at the origin are independent of the coupling constant. On the other hand, ag of course 
changes a lot between the fits in the two different gauges to compensate for the differences 
in the structure of the vector part of the potential, while the other parameters remain almost 
constant. 

The bottom row in Table | shows the total for the four fits. A fictitious error of 10 MeV 
for all masses has been used to calculate the values in order to make them comparable. We 
see that we obtain a slightly better result in the Coulomb gauge and that in the Feynman 
gauge we obtain the largest difference between the fits with different cTj. 

In Table |l|the values in the row just above the lepton pair decays measure the quality 
of the calculation of the El, Ml, and two gamma transitions. We have here excluded the 
decay \I/(2S') rjc'j as it would totally dominate x^- We see that the results in the Coulomb 
gauge changes most and that the Feynman gauge on average gives a better result although 
the Coulomb gauge with cTj = 50 MeV is better than any of the two fits in Feynman gauge. 

In the last row in the same table we find the values for the lepton pair decays. The 
decays of the D states have been excluded in the calculation as they depend on higher order 
corrections not included in our calculations. We see that the Feynman gauge clearly gives a 
better result. 

In conclusion the Coulomb gauge gives better results for the masses, the Feynman gauge 
gives better results for the lepton pair decays and the results for the two gamma. El, and 
Ml transitions are fairly equal, perhaps slightly better in Feynman gauge. The overall 
differences are not so large that we can clearly tell which one of the two gauges gives best 
results. 

The lepton decay widths are normalized such that the value for J/\l/ agrees with the 
experimental value. We see that the other values for the S states become too large. This 
indicates that we lack a fine-tuning of the shape of the potential. A smaller slope of the 
linear confinement potential would push the wave function further away from the origin but 
would on the other hand make the mass spectrum worse. That the decay widths for the D 
states are rather small is not surprising, as they mainly come from higher order corrections 
not included here. In non-relativistic models with a Coulomb potential the wave function 
for D states is zero at the origin, which give a vanishing decay width. In the Salpeter model 
this means that higher order corrections to the kernel in the ladder approximation are very 
important for the D states. 

The values for the two gamma decays agree very well with experiment, and the values 
for the El transitions agree fairly well although the values for the decays including the xa 
state all come out about a factor of two too large. 

When it comes to the Ml transitions we see that the predicted value of the J/"^ —>■ '-frjc 
decay width agrees very well with experimental data. On the other hand the ^^(25) decay 
width is far too large. This is because the wave function of the excited state is not determined 
well enough, and as this decay is forbidden it is extremely sensitive to the exact shape of 
the wave functions. 

For light quarks, the non-relativistic quark model parameters usually relate the magnetic 



12 



moments directly to the constituent masses and there is no room in the parameterization, 
so to speak, to introduce anomalous magnetic moments. 

For the charm quark the current and constituent masses should be essentially the same 
and there is then a possibility to introduce an anomalous magnetic moment as a new inde- 
pendent parameter. 

The details of the decay mechanism are in principle sensitive to this parameter. In 
the absence so far of any measurement of the magnetic moment of charmed baryons, the 
radiative decay rates therefore constitute at present the only place where this parameter can 
be studied. 

Once an accurate determination of the wave functions can be made, the anomalous 
magnetic moment can be calculated from the Ml transitions. But it is important to notice 
that the outcome of such a calculation depends heavily on the quark mass. Therefore 
a determination of the quark mass must be made first. An estimate of the anomalous 
magnetic moment have previously been made in a Salpeter model where a non-relativistic 
approximation has been made p2| . 



A calculation of lepton pair decays and El transitions has earlier been done by Ito 
in a model where only the positive energy term is used in the Salpeter equation. Recently 
a calculation of all decay observables has been done by Resag et al. [^] in a Salpeter model, 
but with a cut-off in the potential to regularize the divergence in the lepton pair decay 
amplitude. The overall results in both these calculations are of about the same quality as in 
ours. As we have discussed earher, the lepton pair decays depend of the cut-off parameter, 
and therefore have to be renormalized before they can be compared. Doing so, we see that 
their results for a scalar confinement are almost identical with our results in Feynman gauge 
with (T, = 10 MeV. 



VIII. IMPORTANCE OF NON-DOMINANT TERMS 



Sometimes the Salpeter equation is reduced to only the term with the positive energy 
states of the quark. This is done by making the following replacement of the propagators 
(0) in the Salpeter equation: 



A+(q)7° 



7°A:(q) 



m 



(gO-E(q)+^e) (g" + E(q) - ze) ^^^^ 
This imposes an extra constraint on the reduced wave function \E'(q). Eq. is replaced by 

Al(q)vl/(q)A;(q) = (40) 
For the wave function only the first term in (10) remains: 

i{M - 2E) 



(M/2 + gO - E + ze)(-M/2 + + E - ie) 



-A+(q)vl/(q)A:(q) 



(41) 



Does this reduction change the results significantly? Olsson et al. |T3[ have recently 
investigated the validity of the reduced Salpeter equation by comparing the eigenvalues and 
wave functions found from the full and reduced Salpeter equation. They found that there is 
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only a significant difference for light quarks, and for as heavy ones as the charm quark the 
difference is very small. 

However, this does not completely test the validity of the reduced Salpeter equation. 
In particular this does not tell us how important the last two terms in (0) are. They do 
not appear in the calculation of the mass spectrum, which only involves the reduced wave 
function \E'(q). On the other hand, the El and Ml transitions do include these two terms. 

To test the validity of the reduced Salpeter equation we have therefore compared the 
decay widths for El and Ml transitions using only the terms with positive energy states 
and the widths using all terms. We have not solved the reduced Salpeter equation but only 
neglected all but the dominant term. This will not give us an exact comparison with the 
reduced Salpeter equation but as have been shown by Olsson et al. [|l^ the wave functions 
differ very little. The result is shown in Table |V[ We see that the decay widths change 
considerably if we neglect all but the dominant term, in some cases with more than 100 %. 
This shows the importance of keeping all terms. Although the magnitude of the negative 
energy part of the reduced wave functions is small - the contribution to the norm is typically 
< 5 % - there are other terms in the decay amplitudes which make the contribution from 
the non-dominant terms significant. 

This can be understood if we note that the last two terms in the wave function (|T^) do 
not appear in in the expression for the norm (|B6|). Using the reduced Salpeter equation 
leads to a normalization condition with only the first term in ( P6|) . The omitted second 
term is small and thus changes the normalization very little. 

To further illustrate this we have plotted the functions 

Tfl 



m 

and 

loiq) 
E{q) 

for the state in Fig. ^. These functions appears in the amplitudes of 0^~'s wave function. 
We see that negative energy amplitude -^^^ p{<i) — ^ A{,<i) is considerably smaller than the 
positive energy amplitude. The amplitude Iq/ E appearing in the last two terms of the wave 
function (|TD|) is somewhere in between in magnitude. However, it should be noted that it is 
difficult to directly compare this function to the other two because of the other factors in 
the amplitudes coming from the propagators. But we see that it is clearly not negligible. 

Comparing the mass spectrum and wave functions is therefore not enough to tell if an 
approximation is valid. The decay widths must also be studied. Using the reduced Salpeter 
equation means that several terms in the decay widths are omitted and we have shown that 
these terms are very important in the El and Ml transitions. We conclude that the reduced 
Salpeter equation gives approximately the same result as the full Salpeter equations only 
for those quantities which do not depend on the last two terms in the wave function (p!oD . 
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IX. SUMMARY AND CONCLUSIONS 



We have done a thorough analysis of the charmonium system in a Salpeter model with 
an instantaneous potential. We have presented general equations for calculating the mass 
spectrum as well as decay widths for lepton pair decays, two gamma decays, Ml and El 
transitions. This has been done for a vector plus scalar potential both in the Coulomb and 
Feynman gauges. The formalism used is based on the one developed by Suttorp |jl|. 

We have discussed the problem with the renormalization of the lepton pair decay am- 
plitude in the instantaneous approximation. This approximation makes the renormalization 
constant Z\ depend on q and Zi — oo as g ^ oo. With a linear confining potential we have 
to make the renormalization in the region where the quarks are asymptotically free and the 
linear potential is negligible. But since Z\ ^ oo this leads to an inconsistency, and therefore 
the only natural way to define the renormalization constant is to normalize the calculated 
widths to the experimental value of the width J/\l/ — >■ e^e~ . 

Numerically we have obtained an excellent fit to the mass spectrum of charmonium. Most 
of the decay widths agree rather well with experiments, although the lepton pair decays of 
the excited S states are too large. The excited states are obviously more sensitive to fine- 
tuning of the model. A smaller slope of the linear confinement potential would push the wave 
function further away from the origin but would on the other hand make the mass spectrum 
worse. Similarly the poor result for the forbidden Ml transition \E'(2S') rjc'j indicates a 
lack of fine-tuning of the potential. We conclude that before any detailed predictions can be 
made from the Salpeter equation, a systematic way to improve the kernel must be developed. 

Comparing the two gauges, Feynman and Coulomb gauge, we have found the the 
Coulomb gauge gives better results for the mass spectrum and that the Feynman gauge 
gives better results for the lepton pair decays. The overall difference is not large enough to 
clearly tell which gauge gives the best results. 

The validity of the reduced Salpeter equation has been tested, and it has been found that 
the decay widths of the El and Ml transitions change considerably if one uses the reduced 
equation. 
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APPENDIX A: WAVE FUNCTIONS FOR DIFFERENT STATES 

We here present the expansion in the Dirac algebra of the reduced wave functions, ^^(q), 
defined in (^, for different states. These expansions are the most general ones with the 
correct transformation properties regarding angular momentum, parity, and charge parity, 
which also fulfill the conditions (^^. All scalar functions \E'j depend only on g = |q|. 

0--^ : ^(q) = (^^p{qh' + ^A{q) (7V + 7^^) ) ^oolq), (Al) 
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+ {^^v{qh + ^TM (a-z7'^)) • Yil(q), (A2) 
: M/(q) = (^p{q)^' + ^A{q) (iV + ^)) Y^U^) , (A3) 



m 

0++ : v[/(q) = (^^^,,{q) + . + vl/^(g)^^ Foo(q), (A4) 

1^+ : ^(g) = (^y(g)^^^ + ^T(g) (.^^ + ^7^7)) ■ Y£(q), (A5) 

2++ : *(q) = (^^s{q) (^1 + ^7 " q) + ^T,(g)^) >^2„.(q) 

+ {^^v{qh + ^TM (a-zf • Ya(q). (A6) 

Here the Y£m(q)'s are the scalar spherical harmonics and the Y£^(q)'s are the electric vector 
spherical harmonics. The latter are defined in terms of the vector spherical harmonics, 

Y7l= E YLm,em,iLm, Im^lJm), (AT) 

mim2 

where e± = ^{e^ ± iey) / \pl and eg = e^, as 



APPENDIX B: NORMALIZATION CONDITION 



The wave functions satisfy the normalization condition 

^ / / fS^^^^^^'^^^^^' ^' + ^' ^^^^^^^ = ^^^^ 



where 



J(p,g,P) = (27r)V^)(p-g)Sri(P/2 + p)52-i(-P/2+p). (B2) 
The adjoint wave function x is defined 

X = 7V+7°- (B3) 



We have contracted with P^^ in (pT|). This is possible if the Salpeter equation is written 
in a covariant form. This is done by letting the kernel depend on the four-vectors (0, q) 
and (0, p) in the rest frame. The Salpeter equation is then Lorentz transformed by first 
multiplying it by the spinor transformation matrices >S'(A) and S'^^(A) from each side. Then 
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the change of variables q h~^q is made. The kernel is found to depend on the orthogonal 
part of the momentum q± = q — j^P and p± = p — in a frame where the particle 
moves with momentum P. In the rest frame q± reduces to (0, q). Thus K = K{q± — p±). 
For such kernels we have 

P'^^K{q^-p^) = 0. (B4) 



Thus the second term in (Bl) does not contribute to the normalization. Integrating the 
other term over p gives the following normalization condition in the rest frame, 

TV j ^x(?)7°x(g)(-y7° + i-m) = 2iM. (B5) 

Using the form ([T0|) of the wave function and integrating over g° gives 

Tr / ^ (vl/+(q)A+(q)vl/(q)A:(q) - v[/+(q)A+(q)vl/(q)A;(q)) = 2M. (B6) 

Note that due to the covariant formulation of the normalization condition, the wave function 
is correctly normalized in any frame. 

By taking the trace in ( p6|) for the different states we find the following normalization 
conditions: 

0-+ : 2M = - r ^q^E{q)^A{q)^piq), (B7) 
m Jo [ZTcy 



1- : 2M = 8 r -^g2 ( :^vl;^(g)vl;^^(g) + :^v[;^^(g)vl/^(g) ] , (B8) 

dq 2; 



1+- : 2M = - —l^^q^Eiq)^Aiq)^piq): (B9) 
m Jo (ztt)'^ 

0++ : 2M = 8 ^gE(g)^5(g)^T(g), (BIO) 

1^^: 2M = -8 rj^^q'^^y{q)^T{q), (Bll) 

2++ : 2M = 8 r -^g^ (^^siq)^TM + ^^TM^viq)] ■ (B12) 
Jo {2ir)'^ \ q m J 



APPENDIX C: EIGENVALUE EQUATIONS FOR THE DIFFERENT STATES 

Equations for the scalar amplitudes of the different states can be found from (|ri|). By 
expanding both sides of (|1T]) in the Dirac algebra, and then multiplying by (Y^mq)* or 
(Y^^(q))* and summing over m, and finally integrating over the angular variables we get the 
equations below [|l|. The functions F/{p, q) and Ff{p, q) are the coefficients of the spherical 
harmonics in the expansion of the kernel given in Section The function F^'^{p, q) is the 
same as F^{p,q), but the terms with F^'^ip^q) only appear in the Coulomb gauge. Thus, 
in the Feynman gauge, these terms should be set to zero. 
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0-+ state 



£;(g)7o (27r) 



state 



M9tM = m f 2 + f 



E(g)yo (27r)3^ V 3p 
1 dp ,(V2mU2 y y jp' - q^) y\ 

+ |^(2Fo^ - Fo^^ + F,^ + F,^^)*^(p)) , (C4) 
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+ 



+ 



1 



E{q) 
1 



dp 
(2^3 

dp 



p 



m 



E{q)Jo (27r) 



2g2 



Fr^nip) - {^Fr + ^(Fo^^ - fD) . (C6) 



1+ state 

m dp ^2 fr.S PI fj^s , oz^5^ 



M^^(g) = m I 2 + 1 ^'p(g) 



0++ state 



1++ state 

M*t(?) ^-171(2 + j 

^ p^(Ff-2F,^)*^(p), (Cll) 



F(g) Jo (27r)3 



M*,(,) = -« + (/^fl-ill)*,te) 
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+ 



m 

W) 

2 



dp 
(2^ 



{2F,' + Fi))^T{p) 



dp 

3E{q)m Jo (2^ 



pq 



2>m? 

V , i^y^ 



T)) *t(p). 



(C12) 



2++ state 

M«T,(«) = m fa + fvl?) 



E{q) Jo (27r)3^ VV 5p ' ' ' ' 2p^ 

2m \ 
+ ^(3(F,^ - F^^) + (2^3^ + 3F3^^))vl>v(p) j , (C14) 



E(q) Jo (27r 
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APPENDIX D: AMPLITUDES FOR THE DIFFERENT STATES 



For calculation of different decay amplitudes we need the expressions for the wave func- 
tions xil) for the different states. Like the reduced wave function, \l/(q), xil) is expanded 
in the Dirac algebra by imposing the correct transformation properties regarding angular 
momentum, parity, and charge parity. The scalar amplitudes of this expansion are then 
expressed in the amplitudes of ^I/(q), by similarly expanding ( [T(]| ) in the Dirac algebra using 
the wave functions given in Appendix 0. 

0"+ state 

The most general wave function for the 0^^ state can be written 

Xiq) = (XA(g)7' + XB{qh\l ■ q) + XciqhV + Xoigh^ct ■ q))roo(q). (Dl) 
The scalar amplitudes Xiil) £^re 

XA{q) = (A+- + A-+)iM/p(g) + (A+- - A-+)^v[;^(g), (D2) 

XB(g) = (A--A++)i/o(g), (D3) 

Xc{q) = -(A+- + A-^)hAq) - (A+- - A-^)^^p{q) 

+(A-- + A++)-^/o(g), (D4) 

XD{q) = (A+- + A-+)-^vl;^(g) + (A+- - A-+)-Lv]/p(g) 

ZTTl Lhi 

+(A- + A++)^/o(g), (D5) 



where 



^ i{M-2E) 

(M/2 + g0-E + ie)(-M/2 + g0 + E-ie)' ^ ' 

A- = + , (D7) 

^++ '[ i'D'R) 

{M/2 + qO - E + ie){-M/2 + q^ - E + ie)' ^ ' 

A"" = - fD9) 

{M/2 + q^ + E -ie){-M/2 + q^ + E -ie)' ^ ' 



The integral /o(q') is 

Io{q) = U^A{q) + \j -^^pHFo'ip,q) + ^Ff (p, q))^ ^ip) 



(27r)3" ^ ^"'"^ q 
/(|p/(F„''(pw)-ife--ff.-))*.(p). (DIO) 
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1 state 

The wave function for the 1 state can be expanded as 

x{q) = {Xi{q) + X2(g)7° + X3(g)(7 • q) + X4(g)(Q; • q))Fi^(q) 

+ (X5(g)7 + X6(g)a + iX7{qh^{a x q) + ixsiqh^icL x 7)) • YS(q). (Dll) 

The amphtudes are 

771 

+ (A++ + A--)— (D12) 
X2(g) = (A++-A--)l/i(g), (D13) 
Xs{q) = (A+- + A-+)^*5(?) + (A+- - A-+)^vi/^,(g) 

-(A++ + A--)-l7i(g), (D14) 
X4(g) = (A+- + A-+)lvi/^,(g) + (A+- - A-+)^*5(g), (D15) 
X5(?) = (A+- + A-+)^^v{q) + (A+- - A-+)^M/T,(g), (D16) 
X6(g) = (A+- + A-+)l*T.(g) + (A+- - A-+)^^y{q) 

+(A++ + A-)-^/2(g), (D17) 
X7(?) = (A++-A--)l72(?), (D18) 
Xs{q) = -(A+" + A-+)^*T.(g) - (A+- - A-^)^^y{q) 

+ (A++ + A--)^72(?), (D19) 

where 

h{q) ^mU^siq) 

dp 



+ 



(27r) 



/2(g) = -t/*r,(g) 
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0++ state 

The wave function for the 0++ state can be expanded as 

Xiq) = iXi{q) + X2(g)7° + X3iq)h ■ q) + X4(g)(a ■ q))>^oo(q). (D22) 
The amphtudes are the same as for 1 but with Ii{q) given by 

1++ State 

The wave function for the 1++ state can be expanded as 

x{.q) = {xiiqh^ + X2(g)7V + xsiqh^-y ■ q) + xMi^cx ■ q))ri„(q) 



+ ix^iqh'l + Xeiqh^cx + iXiiqM x q) + ^X8(g)(q x 7)) ■ Yri(q). (D24) 
The amphtudes are 

Xiiq) = 0, (D25) 

X2(g) = (A---A++)-J^/i(g), (D26) 

X3(g) = (A++ + A-)— 7i(g), (D27) 

mg 

X4(g) = (A---A++)-^/i(g), (D28) 

^g 

X5(g) = -(A+- - A-+)^^y{q) + (A+- + A-+)^^r{q) 

+(A++ + A--)|^/2(g), (D29) 

X6(g) = (A---A++)i72(g), (D30) 

771 1 

X7(g) = -(A+- - A-+)^*y(g) + (A+- + A-+)-^T{q) 

-(A++ + A")-l/2(g), (D31) 
X8(g) = (A+- - A-+)l*v(g) - (A+- + A-+)^*^(g), (D32) 
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where 



- j (f (2< + F^) - 1^ (b^ + e)Fr^ - 2pg<j) vl;^(p). (D34) 



state 



The expansion for the 2^"'" state is the same as for the 1 state with the changes 
Yira —>■ Y2m., ^im- The amphtudes are the same but with Ii{q) and hiq) given by 

hiq) = mU^siq) 



2 I (27,) 



1 1^4^^^''^^^^"^ 



+ (^i^J" + ^ ((/ + g')i^2''^ - '^pqF^")^ ^tAp)) ■ (D36) 



1. Lorentz transformed states 

To calculate the El and Ml transitions we need the wave functions for some states in a 
frame where the particle is moving with momentum P. 

The Salpeter equation in a general frame (j^) is found from the one in the rest frame 
as follows. First this is multiplied by the spinor transformation matrices 5'(A) and S~^{h) 
on each side. A is the Lorentz matrix transforming (M, 0) to the momentum P. Then the 
change of variables q — >■ K~^q is made. 

The Lorentz transformed state is therefore found to be 
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X{q, P) = 5(A)x(A-'?, M)S-\K), (D37) 

where x(g, M) is the wave function in the rest frame. Under the change of variables q — > 
A~^g, then 

q'^qp- ^f, (D38) 

E^uj^ yjq^ + m^, (D40) 

sin ^ ^ sin ^' = -^^^ -, (D41) 

cos ^ ^ cos 9' = ^ , (D42) 

Mqr 

if the particle moves in the z direction. 

For the 0~+ state the transformed wave function is 

X{q, P) = (xAiqh' - XB{qh'i± + Xc{qh'^ + XD{q)^{hP - Ph)^ Foo(^', (D43) 

The scalar amplitudes are the same as before with the above changes of variables, and with 
the functions depending on q-x- 

The Lorentz transformed 1 wave function is 

x(?, P) = {xi + X2^ - X3i± + X4^{4±P - Pi±)j YUO', ^) 

- x.r + §^{rp - PT) + Y^n± - i±r) + xs^{n± - i±n (d44) 

Here 
with 

Y = {Y,{e\ (^), Yy{e\ (^), Y,{e\ ^)). (D46) 

The Lorentz transformed wave functions for 0++ and 2++ are identical to the one for 

1 (only the first four terms for O"^"*"). 
The wave function for 1"'""'" is 

X = (xi7' + X2f ^ - Xaf ^± + XA^{i±P - YimiO', ^) 

- X57T + xe^iTP - m + X7^(r^^ - hf) + xs^ntii. - i±n i^m 
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APPENDIX E: 77 DECAYS 

We here present the matrix element for the 77 decays expressed in the scalar amplitudes 
defined in Appendix ^ We write each scalar amplitude of the wave function as 

Xi{q) = x^(g)A+- + xr^(g)A-+ + x++(g)A++ + xriq)^", (ei) 

where the A^-^ 's are defined in ( P6|) -( p9D . For example in (|D2| ) we write 

XA{q) = xl"(g)A+- + XA+(g)A-+ + xi+(g)A++ + xa"(9)a-, (E2) 

with 

x\- = ^^p(g) + (E3) 

Xa^ = \^p{<D - (E4) 

= Xa' = 0. (E5) 
We then take the trace and perform the integration. This gives the following results. 

0"+ state 

^f^^'^M ( (Xc~(q - k) + xi'mq) ■ (ei X €2) 

eie^iv/v y (27r)3(M-2E(q)-2E(q-k))E(q-k) 

c^q (Xc^(q - k) + xB^mq) ■ (ei x €3) 



+ 16i 
+ 16i 



(27r)3 (M + 2E{q) + 2E(q - k))E(q - k) 
dq 1 



X 



(27r)3 (M - 2E(q) - 2E{q_ - k))(M + 2E(q) + 2E{q - k))E{q - k) 
HXb' -xV)E{ci-k)q 

+ iXc'^ + Xc")(q - k) + + XD")"^q) ■ (d x es)] . (E6) 



16 
16 



0++ state 

rfq (-mxA~ + q • (q - k)xc~)(ei ■ £2) - 2xc~(ei • q)(e2 • q) 
(27r)3 (M-2E(q) -2E(q-k))E(q-k) 

dq. {-mxA^ + q ■ (q - k)xc^)(ei ■ €2) - 2xc^(ei ■ q)(e2 ■ q) 



(27r)3 
dq 



(M + 2E(q) + 2E(q - k))^(q - k) 
1 



X 



(27r)3 (M - 2E(q) - 2E(q - k))(M + 2E(q) + 2E(q - k))E(q - k) 

HXa^ + Xa") + ^(q - ^){Xb' - Xt^) 



+ q- (q-k)(xJ+ + Xc7 ))(ei-e2)-2(xJ+ + Xc ) (^i " q) (^2 ■ q) 



(E7) 
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I 



dq 



2++ state 
1 



(27r)3 (M - 2£;(q) - 2E{q - \s.))E{q - k) 



X 



{-mxt +q-(q-k)x3 )(ei-e2)-2x3 (ei-q)(e2-q) 

+ Xt~ ■ (q - k)(ei • €2) - (ei • xt~)i^2 ■ q) - (ei • q)(e2 • xt 
dq 1 
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r _dq 
J (27r) 



X 



(27r)3 (M + 2E{q) + 2E{q - \s.))E{q - k) 
(-mxr^ + q • (q - k)xj+)(ei • es) - 2^3 +(ei • q)(e2 • q) 

+ X5^ ■ (q - k)(ei • €2) - (ei • X5^)i^2 ■ q) - (ei • q)(€2 • X5 
dq 1 
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X 



(27r)3 (M - 2E{q) - 2E{q - k))(M + 2E{q) + 2E{q - k))E{q - k) 
{-mixf + XD + E{q - k)(x2" - xt^) 

+ q • (q - k) {xt^ + X3"")) (ei • 62) - 2(xJ+ + X3"") (^i • q) (^2 • q) 
+ (xr + Xr) ■ (q - k)(ei • €2) - ((61 • + X5"))(e2 • q) 



+ (£2 • + ))(ei-q)) 



(E8) 



Here Xbiq) = X5(9)YS(q)- 
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FIGURES 



FIG. 1. Feynman diagram for lepton pair decays. 
FIG. 2. Feynman diagrams for two gamma decays. 
FIG. 3. Feynman diagrams for electromagnetic transitions. 

FIG. 4. Amplitudes for the ^ state. We see that the function — appearing the 

negative energy part of the amphtude is much smaller than the function ^^p + in the positive 
energy part of the amplitude. However, the function Iq/E appearing in the mixed energy part is 
somewhere in between and certainly not negligible compared to the positive energy amplitude. 
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TABLES 



TABLE L Masses of the charmonium particles. The four columns in the middle contain the 
values obtained in the four different fits. F stands for Feynman gauge and C for Coulomb gauge, 
10 MeV and 50 MeV stands for the minimum value of the experimental error, cjj used in the fit. 
We see that all four fits give excellent agreement with experiments. In the last row we give 
including all masses using an fictitious error of 10 MeV for all masses in order to make the numbers 
comparable. We see that the Coulomb gauge gives a slightly better result. Experimental values 



are taken from Review of Particle Properties data table |24]. 



particle 


F 10 MeV 


F 50 MeV 


C 10 MeV 


C 50 MeV 


exp 




2.97 


2.97 


2.97 


2.97 


2.98 


r/c(25)[not used in fit] 


3.62 


3.64 


3.62 


3.63 


3.59 


J/^ 


3.14 


3.15 


3.13 


3.14 


3.10 


*(25) 


3.70 


3.72 


3.69 


3.70 


3.69 




3.76 


3.77 


3.76 


3.77 


3.77 




4.09 


4.10 


4.09 


4.10 


4.04 


^(3D) 


4.13 


4.13 


4.14 


4.14 


4.16 


^(45) 


4.39 


4.39 


4.40 


4.39 


4.42 


XcO 


3.44 


3.41 


3.43 


3.41 


3.42 


Xcl 


3.50 


3.51 


3.50 


3.52 


3.51 


Xc2 


3.49 


3.50 


3.49 


3.51 


3.56 


hci [not used in fit] 


3.49 


3.50 


3.49 


3.51 


3.53 


x' 


127 


153 


112 


116 
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TABLE II. Predicted values for the decay widths compared to the experimental values. The 
width J/^' — > e^e~ is used as input for the renormalization. The row with values above the 



lepton decays presents the total for the two gamma, El, and Ml transitions, except ^'(25') 
which would otherwise totally dominate x^- The last row is the total for the lepton decay, but 
with the decays of the D states excluded, as they depend on higher order corrections not included 
in our calculations. Experimental values are taken from Review of Particle Properties data table 
unless otherwise stated. 



decay 




F 10 MeV 


F 50 MeV 


C 10 MeV 


C 50 MeV 


exp 


r/c 77 




6.2 


6.3 


6.2 


6.5 


7.5 ± 1.5 


XcO 11 




1.6 


1.8 


1.5 


1.6 


1.7 ± 1.3^ 


Xc2 11 




0.31 


0.41 


0.31 


0.37 


0.37 ± 0.17 


J/^ r/e7 


1.65 


1.33 


1.66 


1.38 


1.14 lb 0.39 


^{2S) 


Vcl 


9.8 


10.9 


10.2 


12.8 


0.78 ± 0.24 


^(25) ^ 


V'cl 


0.27 


0.17 


0.26 


0.18 




XcO ^ J/^1 


130 


96 


143 


110 


92 ± 40 


Xcl J 1^1 


390 


399 




AO A 

434 


240 ± 40 


Xc2 J/^1 


218 


195 


240 


218 


267 ± 33 


■^{2S) 


XcOl 


31 


47 


26 


31 


26 ±4 


^{2S) 


Xcll 


58 


49 


63 


50 


24 ±4 


■^{28) 


Xc2l 


48 


47 


51 


49 


22 ±4 






136 


127 


174 


116 




J/-^ e 


^e- 


5.26 


5.26 


5.26 


5.26 


5.26 ±0.37 


^{2S) 


e+e" 


2.81 


2.45 


2.92 


2.68 


2.14 ±0.21 


^{2D) 




0.09 


0.19 


0.25 


0.35 


0.26 ±0.04 


^(35) ^ 


e+e~ 


1.99 


1.59 


2.09 


1.84 


0.75 ±0.15 


^{2,D) 


e~^e~ 


0.14 


0.29 


0.05 


0.07 


0.77 ±0.23 


^(45) ^ 


e+e^ 


1.42 


1.00 


1.56 


1.31 


0.47 ±0.10 


x' 




169 


62 


212 


130 





^From Ref. 125 
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TABLE III. Parameter values in the four fits and the value of as defined in Eq. (38). When 
changing o", from 10 MeV to 50 MeV the parameters which change most are the quark mass m 
and the constant U. The parameter naturally changes significantly between the fits in the two 
different gauges, while the other parameters remain almost constant. We see that is slightly 
smaller in the Coulomb gauge when cij = 10 MeV and slightly smaller in the Feynman gauge when 
a,; = 50 MeV. 



Parameter 


F 10 MeV 
= 110 


F 50 MeV 
= 10.2 


C 10 MeV 

X' = 95 


C 50 MeV 
= 11.3 


as 


0.282 


0.288 


0.348 


0.369 


m 


1.366 


1.276 


1.358 


1.298 


A 


0.257 


0.280 


0.266 


0.274 


U 


-0.131 


0.072 


-0.112 


0.070 



TABLE IV. Change of decay amplitudes when only the term with positive energy states is 
used. We see that most decay widths would increase, some with several 100%, when only the term 
with positive energy states is used. This shows that it is important to keep all terms in the model. 



Decay 


F 10 MeV 


F 50 MeV 


C 10 MeV 


C 50 MeV 




38% 


66% 


40% 


60% 




-47% 


-67% 


-59% 


-68% 


^{2S)^ri',j 


117% 


240% 


117% 


212% 


Xco J/*7 


185% 


357% 


176% 


307% 


Xci J/*7 


15% 


32% 


10% 


21% 


Xc2 ^ 


58% 


119% 


54% 


97% 


^{2S) ^ Xco7 


65% 


100% 


68% 


101% 


^{2S) ^ Xci7 


-17% 


-24% 


-15% 


0% 


^{2S) ^ Xc27 


21% 


33% 


18% 


27% 
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